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ABSTRACT 

The core of planetary nebula NGC 6302 is filled with high-excitation photoionised gas at 
low expansion velocities. It represents a unique astrophysical situation in which to search 
for hyperfine structure (HFS) in coronal emission lines from highly ionised species. HFS 
is otherwise blended by thermal or velocity broadening. Spectra containing [A1VI] 3.66pm 
3 ?2 <— 3 Pi, obtained with Phoenix on Gemini-South at resolving powers of up to 75000, 
resolve the line into five hyperfine components separated by 20 to 60 km s -1 due to the 
coupling of the 7 = 5/2 nuclear spin of 2 ‘ Al with the total electronic angular momentum J. 
26 A1 has a different nuclear spin of I = 5, and a different HFS, which allows us to place a 3 cr 
upper limit on the 26 Al/ 2 ' A1 abundance ratio of 1/33. We measure the HFS magnetic-dipole 
coupling constants for [Al vi], and provide the first estimates of the electric-quadrupole HFS 
coupling constants obtained through astronomical observations of an atomic transition. 

Key words: atomic data - atomic processes - linefidentification - line: profiles - ISM: abun¬ 
dances - platenary nebulae: NGC 6302. 


1 INTRODUCTION 

The interaction between the electronic wave-function and a non¬ 
zero nuclear magnetic dipole splits a fine-structure level { L , J} 
into hyperfine levels. While hyperfine transitions are common in 
the radio range, at shorter wavelengths atomic hyperfine struc¬ 
ture (HFS) has seldom been resolved in emission. Examples of 
hyperfine broadening include solar absorption lines from neutral 
species ( Abll 19521 ). such as Mn I, whose transiti ons at 1.7743 jam 
are a rare example o f unblended hyperfine lines l Melendez 1 19991) . 
[Booth & Blackwell (and references therein, Il983) ~ summarise the 
effects of HFS on stellar spectra: its neglect results in incorrect 
measures of line broadening, and mismatched curves of growths, 
leading to ~0.2 dex errors in the inferred photosphere elemental 
abundances, even for faint lines far from saturation. Both HFS and 
saturation lead to line broadening. This degeneracy complicates the 
use of stellar absorption lines as diagnostic of the hyperfine cou¬ 
pling constants. Hyperfine splitting can also be resolved in inter¬ 
stellar Na I D-line absorpti on in the local ISM , which requi res re- 
solving powers of > 5 10 5 IWavte et alll978l : |Barlow et alll995l 
for more recent data). 

The 1.8 MeV gamma-ray emission due to the decay of 26 Al 
into 26 Mg has been the object of extensive space borne surveys: 
with a half-life of 7.2 10 5 yr, 26 Al is a signpost of recent nucle¬ 
osynthesis. Line emission at 1.8 MeV from the diffuse ISM (as 
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observed by COMPTEL. iDiehi e t al|[l995i) is consistent with an 
26 Al source in either AGB stars"^^^rtin^t al| l99ll) . novae, su¬ 
pernovae, Wolf-Rayet stars l Pr antzo si [2004 ) or from cosmic-ray 
c ollisi ons jn molecular clouds dClavtonlll99i) . The INTEGRAL 
fWinkler et all2003l) mission holds the promise of improved angu¬ 
lar resolution with which to identify the most important contributor 
to the diffuse emission. 

Although the decay of 26 Al is observed in the ISM at large, 
the 26 Al / 21 Al isotopic abundance ratio (hereafter Ri so ) has never 
been measured in any astrophysical source. The only available up¬ 
per limit in any specific object is that of lBaneriee et all 12004) . who 
observed the vibronic bands of AlO at 1.5 /im in the nova-like vari¬ 
able V4332 Sgr, and reported an upper limit of ~1/10, lacking a 
statistical discussion. 

As an application of our detection of HFS in [Al VI] 3.66/rm 
3 P 2 <— 3 Pi (hereafter [Al VI]), which is the first in an astrophysical 
near-IR emission line, we can set an upper limit on 26 A1/ 27 A1 using 
the difference in the HFS of both isotopes: the stable isotope 2 'A1 
has a nuclear spin I = 5/2, while 26 Al has a nuclear spin 7 = 5. 
This new upper limit is the most stringent obtained so far in any 
astrophysical target. 

The first detection of atomic HFS in emission, aside from 
the 21 cm Hi line, is to our knowledge the observation of re¬ 
solved HFS in r 13 C III 157.8n m 2 P 1/2 +- 2 P 3/2 by [Boreiko et alj 
1 1988). iKellv & Lacvl i 1 9951) identified multiple components in 
[Na IV] 9.0pm with the hyperfine splitting of 3 P 2 <— 3 Pi. Although 
[Naiv] 9.04/im and [Alvi] 3.66pm are the same fine-structure 
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transitions from isoelectronic ions, they differ in nuclear spin and 
electronic wave-functions, leading to different hyperfine structures. 

We also derive values for the electric quadrupole constants 
(hereafter B constants) in the [Alvi] transition. To our knowl¬ 
edge, this is the first measurement of such constants in an 
atomic transition in any astrophysical object, although the B 
quadrupole constants have been measured in molecular transitions. 
In contrast with atomic HFS, in molecules the hyperfine split¬ 
ting of a given rotational transition primarily derives from nuclear 
quadrupole moments rather than from nuclear magnetic moments 
(e.g.. lTownes & Schallowl 195j). For instance B val ues have pre¬ 
viously been measure d by Tumer & Gammo3 1 1 9751 ) in CN(K = 1 
-0) at 2.6 mm, and bv lziurvs et alJ <1992h in HCNH+(J = 1-0) at 
74 GHz. In this work we show that the inclusion of the atomic elec¬ 
tric quadrupole terms has important spectroscopic consequences. 
It allows improved measurement of the magnetic dipole coupling 
constants by lifting the statistical bias between the magnetic dipole 
constants of the upper and lower levels. 

In this work we demonstrate the use of HFS itself as a di¬ 
agnostic tool in the context of planetary nebulae (PNe). Atomic 
hyperfine effects have previously been used bv lClegg~ akl <1997t) 
in ClII] A1909.6 <- 3 Po to measure 13 C/ 12 C in PNe. 

They recognised that the non-zero nuclear spin of 13 C addition¬ 
ally 1 mixes the 3 Po and 3 Pi fine-structure states. Cm] A1909.6 
is dipolar-electric in 13 C, while it is completely forbidden in 12 C 
because it has no net nuclear spin. This C III multiplet is thus com¬ 
posed of three lines, one of which is due solely to 13 C. 

NGC 6302 is the highest excitation PN known, with a spec¬ 
trum rich in molecular lines, dust, and coronal ions such as 
[Si IX] 3.93 /rm, which can only be produced by photons harder 
than 303 eV, or by electron collisions at T e w 10 6 K. Its 
spectrum can be reproduced by ionisation-bounded photoio nisa- 
tion models with a T = 250 000 K central star JCasassus et alJ 
1200(1) . and the absence of a fast wind makes improbable a signif- 
icant c ontribution from shock excitation. Although the report of 
iMeabum & Walshlll98ct ) for broad wings under [Ne V] 3426 Ahas 
been taken as evidence for a fast wind in NGC 6302, a 3000- 
10000 A echellogramme we acquired with UVES on the VLT 
(C asassus et al., in pre paration) does not confirm the observations 
of IMeabum & Walshl Il98(ll . The photoionised coronal li nes in 
NGC 6302 are astonishingly narrow lAshlev & HvlanJl988l l com¬ 
pared to conditions of collisional ionisation where their abundance 
is maximum. The line-widths measured bv ICasassus et all t 2000 ) 
reflect negligible thermal broadening from photo-ionised gas tem¬ 
peratures of 20000 K, and very small expansion velocities in a 
filled-in nebula. 

It is its small expansion velocity and rich spectrum that makes 
NGC 6302 an ideal object for the use of hyperfine structure as a 
diagnostic tool. 

We describe data acquisition in Section [2] then data analysis 
and results in Section [3] and summarise our conclusions in Sec- 
tion[4] Data reduction and analysis were carried-out using the Perl 
Data Language l http : //pdl. perl. orgl. 


2 OBSERVATIONS 

We observed NGC 6302 with Phoenix fainkle et al J2003j) on Gem¬ 
ini South on 5 nights of May and July 2003, as summarised in 

1 This multiplet arrises from the mixing of 3 P and 1 P states due to mag¬ 
netic interactions between the electrons 


V 

• > 

Figure 1 . Overlay of the Phoenix 14 " slit on the R-band Gemini acquisi¬ 
tion image. 

Table Q The slit position angle was 70 deg East of North, and it 
was centred on NGC 6302’s radio core at J2000 RA=17_H3:44.4, 
DEC=-37:06:11.2, as inferred from the 5 GHz map of lGomez et al J 
i 19891) . at the position of the intensity decrease in the centre of the 
putative radio torus. Fig. Q shows the slit position overlaid on the 
R-band image obtained with Gemini's acquisition camera. Back¬ 
ground cancellation was obtained by differentiation with a refer¬ 
ence field devoid of nebular emission, offset 40" North of the neb¬ 
ular core. Typical integration times in the [Al VI] settings were lh- 
2h on-source for each night, but the noise level largely reflects the 
weather conditions. Poor weather also results in inaccurate back¬ 
ground cancellation. The seeing has a direct impact on the reso¬ 
lution of the spectra, by convolving the emission in the slit with 
neighbouring emission from the expanding nebula. The emission 
that falls through the slit is the convolution of the slit aperture with 
the point spread function. Because of the spatial variations of radial 
velocity within the nebula, poor seeing allows emission from mate¬ 
rial with a wider range of velocities to be admitted by the spectro¬ 
graph slit. The resulting spectra are therefore degraded by a com¬ 
bination of the instrumental resolution and the spatial variations in 
velocity. 

The acquisition of a precise position in the nebula is impor¬ 
tant to obtain consistent spectra. To centre the slit on the position 
of NGC 6302’s radio core we peaked-up on a reference astrometric 
standard in the K band, then offset to the object, switched detector 
settings and applied an additional offset to account for the differ¬ 
ence in refraction between the filters. Peaking-up with a narrow 
slit is difficult because of variations in the seeing on time-scales 
comparable to the acquisition procedure. The accuracy involved 
in peaking-up depends on the seeing and the slit-width. The over¬ 
all positional uncertainty is ~0.6 arcsec, as estimated by adding 
in quadrature the errors involved in peaking up, of about 0.35" 
(or twice the slit-width), in offsetting from the reference star, of 
^0.5 ", and in the filter change, of ~0.l" (Bernadette Rodgers, 
private communication). 

The resulting spectra for the three best nights are shown in 
Figja after subtraction of a small level of continuum nebular emis¬ 
sion. Wavelengths are given in air and in the observatory rest frame. 
There are at least four e mission features observed near the [Al VI] 
wavelength given bv ICasassus et alJ < 12000 3 of 3.659 /jm. The three 
brightest features share simular position-velocity structures. 

Because the 1 arcsec uncertainty in acquisition is larger than 
the slit width, spectra taken on different nights may not sample the 
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same region in the nebula, so we only coadded the frames taken 
with the same instrumental setting (i.e. within the same night). An¬ 
other reason to avoid averaging all the spectra is to keep a good 
spectral resolution: it can be inferred by inspection of Fig. [2] that 
the night of July 30 th has the best line contrast, even though we 
used the widest slit. The sharpest lines should be obtained with the 
narrowest slit. The quality of the spectrum from July 30 th reflects 
that it was acquired under the best weather conditions. 

No reference lines were visible in the calibration exposures 
taken with a Th-Ar-Ne arc lamp in the [Al VI] instrument config¬ 
uration. Instead we used emission lines from coadded sky spectra 
extracted from the science observations (without differencing the 
nodded frames). We fit a straight line to the position of sky emis¬ 
sion features present in a model high-resolution sky spectrum based 
on the HITRAN database (Rothman et alJI 1992h . The accuracy of 
the inferred dispersion law is checked a-posteriori by comparing 
different nights and previous wavelength measurements. The raw 
spectra are modulated by the atmospheric transmission (hereafter 
AT) spectrum, which is reasonably smooth near [Al VI] (AT does 
not show deep troughs). We correct for the AT modulation by divid¬ 
ing the object frames by the spectrum of a standard star (HR 6789) 
grown along the slit. 

The optimal aperture for spectrum extraction in the spatial di¬ 
rection along the slit was determined by varying the upper y up and 
lower yi 0 rows of detector pixels. We summed all the signal in the 
detector within the rows y\ 0 and y up , and estimate the noise level 
a-posteriori, from the rms dispersion of the output spectrum in a re¬ 
gion devoid of line emission. A search in the 2-D parameter space 
{yio, 2/up} for the best signal-to-noise spectrum gives the optimal 
aperture indicated in Fig. [2] We hereafter refer to spectra extracted 
with this optimal aperture as collapsed spectra. 


3 HYPERFINE SPECTRAL FITS 


In Russell-Saunders coupling the magnetic field due to the nuclear 
spin splits a given { L , J} fine-structure level into hyperfine levels, 
with the following energy shifts iplass & Hibberil97ij) . 


AE(L,J,F,I) = 


+hBL,j 


K(K + 1)--I(I + 1)J(J + 1) 


( 1 ) 


where 


K = F (F + 1) — I (I + 1) — J (J + 1), (2) 

where L stands for the electronic orbital angular momentum 
and I is the nuclear spin, h is the Planck constant. Al,j and 
Bl,j are the magnetic-dipole and electric-quadrupole hyperfine 
coupling constants, respectively. F (F + 1) is an eigenvalue 
of F 2 , where F is the vectorial I + J operator. The relative 
intensities 5({Ji, Fi}, {J 2 , F 2 }) of each hyperfine component 
{I, Ji,F\} <— {/, J 2 , F 2 } can be derived from 

S({I,J 1 ,F 1 },{I,J 2 ,F 2 }) = 

{2F , +1x2a+ 1 ) | j 2 ;}’■ < 3 > 

with the selection rule: 

\Fi — F 2 I ^ 1 ^ Fi + F 2 , (4) 

where { } is the six-j symbol defined bv lBrink & Satchlerl(ll994) . 



Figure 2. Phoenix detector array after flat-fielding, correcting for the slit 
tilt, and background and nebular continuum removal. Intensity is given in 
units of the noise on a linear grey scale covering the full range of intensities, 
x-axis is wavelength in microns, y-axis is offset along the slit in arcsec. The 
y-axis increases towards the East. The horizontal arrows limit the optimal 
extraction aperture. 


The model hyperfine structure of [Al VI] given in Table ^de¬ 
rives from an ab-initio calculation of the HFS coupling constants 
(carried out by o ne of us, P JS) using the ato mic structur e code SU- 
PERSTRUCTURE lEissner et alll974lciegg et all 19971) : 

A% i( 27 Al) = 0.2, A% 2 ( 27 Al) = 3461.8 , 

A t J h =1 ( 26 Al) = 0, A t J h =2 ( 26 Al) = 1333, w 

where all values are given in MHz. Table 0 does not include 
the electric-quadrupole terms in the hyperfine energy shifts, be¬ 
cause SUPERSTRUCTURE does not currently predict the electric- 
quadrupole hyperfine coupling constants. The velocities in Table[2] 
are all given relative to the strongest component of the stable iso¬ 
tope. The nuclear magnetic dipole moments used in the calcula¬ 
tion are + 3.64151 un fo r 27 A1 lRaghavanll 989). and +2.804 pn 
for 26 Al l lCooper et alJll99d) . in units of the nuclear magneton 
Pn = ehl2m v c. 

The hyperfine coupling constants given above were calculated 
in a simple two configuration atomic model, 2s 2 2p 4 and 2p 6 . Cal¬ 
culations were also made for more elaborate configuration bases 
and with different orbital optimisation procedures leading to a 
range of values for the hyperfine coupling constants. The results 
were all within 100 MHz of the values quoted above but with no ob¬ 
vious convergence to one particular best result. We therefore quote 
the results of the simplest calculation and adopt cr th = 100 MHz as 
the likely uncertainty in the theoretical result. This uncertainty has 
to be compared to the difference IAj 1 ^ — Af =2 \ ~ 1000, which is 
roughly how these quantities enter the expression for the hyperfine 


© 2002 RAS, MNRAS OOP. miTl 































4 S. Casassus et al. 


Table 1 . Observation log. S/N is the ratio of the peak specific intensity to the noise in the image. 


date 

2003 

slit 

arcsec 

S/N 

integration 

s 

airmass 

seeing 1 

weather 1 

May 7 

0.25 

14 

84 x 180 

1.27 -> 1.03,1.05 -r 1.4 

2 

2 

Jul 26 

0.17 

12 

40 x 300 

1.10 -* 1.01 -r 1.17 

3 

2 

Jul 27 

0.34 

10 

12 x 300 

1.45 -* 2.08 

3 

3 

Jul 30 

0.34 

12 

12 x 300 

1.28 -> 1.64 

1 

1 

Jul 31 

0.34 

14 

24 x 300 

1.08 -> 1.01 -r 1.02 

2 

4 


1 refers to a relative quality number, assigned by inspection, in which 1 is best. 


Table 2. SUPERSTRUCTURE ab-inido calculation of velocity shifts relative 
to the strongest component of 27 A1, and corresponding relative intensities, 
for [ 26 Al VI] and [ 27 A1VI]. 


27 Al VI 


rel. 

rel. 

vel. 

int. 

(kms- 1 ) 

0.0 

5.0(1)! 

-58.0 

4.0 (2) 

-103.1 

3.0 (3) 

-135.2 

2.0 (2) 

-154.5 

1.0(1) 


26 Al VI 


rel. 

rel. 

vel. 

int. 

(km s _1 ) 


-14.8 

9.0(1) 

-49.5 

7.8 (2) 

-79.3 

6.6 (3) 

-104.1 

5.4 (2) 

-107.3 

4.2(1) 


1 the number of sub-components blended together in each velocity compo¬ 
nent. 


energy shifts. The theoretical A th agree with observations within 
2<r th (see below). The high accuracy of the calculation for 27 A1 
should be carried over to 26 Al since they have the same electronic 
wave functions (to a very good approximation). Therefore it is rea¬ 
sonable, within a < 10% uncertainty on the HFS coupling con¬ 
stants, to use the theoretical hyperfme constants for 26 Al in a fitting 
procedure to look for evidence of 26 Al. 

There is no allowance for isotopic mass shift neither in the ab- 
initio calculations nor in the spectral fits. We assume both 26 Al and 
27 A1 share the same line centroid. The normal mass shift (NMS) 
due to the difference in Rydberg constants between 26 Al and 2 'Al 
would cause the centroid of the 3.6/rm line in 26 Al to be shifted 
to the red by 0.23 km s -1 relative to 2 'Al. The specific mass shift 
(SMS) is not known for [Al VI] but measurements have been m ade 
for the same transition in the isoelectronic Ol iDe Natale et alJ 
[ 1995), which show that the ratio of the total isotope shift to the 
normal mass shift ((NMS+SMS)/NMS) is 1.40 between 17 0 and 

le O and 1.26 between 18 0 and 17 0. Adopting the larger of the 
these two values, we can estimate that the total isotope shift of the 
26 Al centroid relative to 2 'Al should be no more than 0.32 km s _1 . 

For both isotopes there are actually nine lines which in prac¬ 
tice reduce to five due to degeneracy in the hyperfme levels asso¬ 
ciated with the J = 1 state. This degeneracy is not exact but in 
practice the lines lie within less than 0.1 km s -1 of each other, 
which is much smaller than the typical linewidth of cr ~ 8 km s _1 , 
so they can be taken to have the same velocity shift. The number 
of components within each line is given in brackets in Table[2] We 
stress that Table[2]does not include the electric quadrupole terms in 
the hyperfme energy shifts. 


We fit the [Al VI] line profile F\ with the following 
parametrised model, 


-^isotope ^gauss 

F\ = F 0 + E E E 

i= 1 g=l F 1 ,F 2 


Ri R g exp —- 


1(A-A 9 (F 1; F 2 )) S 


( 6 ) 


where F a is a constant baseline, Ni so t op e is the number of isotopes 
(i.e. 1 or 2 in this case), JVg auss is the number of Gaussians used to 
represent the fit (1 or 2), Ri is an overall amplitude for isotope i, R g 
is the relative amplitude of additional Gaussians relative to the first. 
Thus the first Gaussian component for isotope i has R g =i = 1, and 
Ri is constant for all g components of isotope i. X g (Fi, F 2 ) is the 
Gaussian centroid of each hyperfme component, 


A 9 (Fi,F 2 )= (7) 

c h 

Ac + AA 9 + AE ^ Jaj p 2jI ^_ A£(Lj JliFlt iy (8) 

where A 0 is a reference wavelength which does not necessarily 
match the fine-structure transition, which instead corresponds to 
the overall centroid of the line, i.e. the average of each HFS compo¬ 
nent weighted by its flux. AA 9 is an offset to describe the velocity 
profile with Gaussian g (one Gaussian has AA = 0). 

The optimisation was carried-out by minimising \ 2 * * = 

y (F(Xj) — F m (A_,)) 2 /erf- in two steps. We perform an ini¬ 
tial heuristic search of the global minimum with the pikaia ge¬ 
netic algorithm of ICharbonneaul C l995t) . and then optimise with 
the variable-metric routine MIGRAD of the Minuit package 
from ICERNl ( 19981) . We took the precaution of cross-checking 
the MIGRAD results with the downhill simplex method amoeba 
fe*ress et alJll986 ). which we observe to be much slower and far 
less robust than MIGRAD (amoeba requires fine tuning of the in¬ 
put simplex and tolerance parameters). Errors on individual param¬ 
eters are estimated by searching parameter space for the A \ 2 = 1 
contour. 

The resulting observed spectrum and model line profile are 
shown on Fig. [3] We show the case of the optimal-extraction spec¬ 
trum from July 30 th in Fig. ^Ji. together with an indication of the 
hyperfme splitting of 26 Al, had it been present. The coadded spec¬ 
trum is compared to the fits in Fig. 05. where it can be appreciated 
that the inclusion of the electric quadrupole hyperfme terms im¬ 
proves the fit. It can be verified by inspection that the solid line, 
with B terms, is appreciably closer to the data than the dotted line, 
without B terms. The residuals are shown on Fig. with the for¬ 
mal 26 Al fit. 

Our best-fit line profiles are summarised in Table0 The 20 Al 
hyperfme coupling constants were kept fixed at their theoretical 
values, as specified in Eq.0 We list reduced \ 2 as an indicator of 
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goodness of fit: values much less than 1 reveal that we are fitting the 
noise with an excessive number of free parameters. We nonetheless 
include these fits in the list with the goal of combining the results 
from all nights. 

We can use the information that the HFS coupling constants 
are the same on each night to perform a second run of the fitting 
procedure, and fix the HFS constants to the average given in Ta¬ 
ble 0 This allows improved estimates of the fine structure cen¬ 
troid, as well as tighter limits on the abundance of 26 Al relative 
to 2 'A1. The results of this second run of fits are summarised in 
Table0 The average value for the fine structure centroid includes 
a correction for the heliocentric syste mic velocity of NGC 6302 
of —35.0 km s” 1 ICasassus et al jEoOtil) . and is accurate to within 
1 km s" 1 . The constraints we can place on the Al isotope ratio are 
summarised under columns f?f so and 7?f so . The 1 a values and up¬ 
per limits indicated under 7?,“ 0 are formal indicators of the relative 
limits set by each spectra, and are derived from the first fitting pro¬ 
cedure, with free HFS constants. In this case we cannot use the error 
estimates derived from the A^ 2 = 1 contour because the positivity 
requisite on Ri so precludes reaching the global \ 2 minimum with 
certainty. But we relaxed the positivity constraint in the second fit, 
fixing the HFS constants. The results are listed in Table 0 under 
the column The weighted average for 100 x Ri so is 0.6±0.8. 

We estimate an upper limit on Ri BO by generating a synthetic 
spectrum with Ri so = 0.03, and repeating the fitting procedure, for 
100 different realisations of Gaussian noise, at the same level as that 
of the collapsed spectrum for 30-06-2003 (which has the best S/N). 
This Monte-Carlo error analysis shows we can recover the input 
isotope ratio at 2 -a: 100 x Ri so = 3.1 ± 1.6. Another simulation 
with the noise level of the residual spectrum shown on Fig.0gives 
100 x Ri ao = 3.0 ± 0.8. These simulations and the combined 
measurement of Ri BO from Tableware in agreement, which allows 
us to place the following 3-a upper limit: 

Riso < (Riso) + 3a = 3.0 1 (T 2 . (9) 

The uncertainty in the measured HFS coupling constants, rel¬ 
ative to the optimal value, is rather large compared to that of the 
central wavelengths. This is due to a significant statistical bias 
in the values of Aj= 2 ( 2 'A 1 ) and Aj=i( 2 'A1). Fig. 0 is a 2-D 
slice in parameter space showing the correlation of both constants. 
With the neglect of the quadrupole HFS constants Bl,j, the bias 
is much stronger and the uncertainty on the magnetic dipole con¬ 
stants is much larger (~ 5 times larger): Aj = i( 27 A1) = I 49 I 444 , 
Aj = 2 ( 2 ' Al) = 3499j)g 24 . In the absence of the electric quadrupole 
terms, the hyperfine energy shifts depend on the Al, j constants ap¬ 
proximately through their difference, Aj = 2( 27 A1) — Aj = i( 27 A1). 
Notwithstanding this difficulty, the observed constants are close to 
the theoretical values used to produce Table0 

The reasons why we are confident on our detection of the 
electric-quadrupole hyperfine splitting are as follows. 

(i) The fit to the line profile significantly improves, with reduced 
X 2 increasing by more than 0.1 without the Bj( 27 Al) constants, 
systematically for all nights. For example, in the case of the col¬ 
lapsed spectrum for July 30 th , \ 2 l (i) (ii) (iii) (iv) v rises from 1.08 to 1.20. 

(ii) The improvement in the fit is visible to the eye by inspection 
of Fig. 0 where it can be appreciated that the fit without Bj( 27 Al) 
gives an excess in the blue. 

(iii) The two electric quadrupole constants are measured with 
accuracies of 6a and 10a. 

(iv) The average isotope ratio without the quadrupole terms is 

(—1.41 ± 0.59) x 10 -2 , which reflects a tendency to compensate 



Figure 3. a) Points: collapsed spectrum of [Al VI] from 30-06-2003, with 
the optimal extraction aperture. Solid line: the best fit with two Gaussians 
per component, and the parameters given in Table[3] without contribution 
from 26 Al. Grey solid line: the profile of 26 Al. had it been present at a 
level giving an isotope ratio of 1. b) Points: coadded spectrum. Solid-line: 
combined model. Dotted-line: combined model without electric-quadrupole 
hyperfine splitting, c) Histogram: binned residuals, excluding the 26 Al fits. 
Solid line: combined 26 Al fit. 

for the misfit with a negative, and spurious, amplitude for the rare 
isotope. 

The intrinsic profile of [Alvi] is manifestly very narrow. A 
single Gaussian fit to the spectrum from 30-06-2003, with a 6- 
row extraction centred on the peak of emission along the slit, 
gives a width of a = 7.31 ± 0.24 km s _1 , or a FWHM of 
17.21±0.56 km s , This spectrum was acquired with the widest 
slit, and a resolving power R = 40 000. We cannot give a precise 
measure of the instrumental resolution because of the lack of arc 
lines near [Alvi], But assuming 3 10 5 /R km s _1 corresponds to 
the FWHM instrumental resolution within 20%, we can give an es¬ 
timate of the de-convolved line width of 15.5±1.0 km s -1 FWHM. 
A comparison with emission lines from lighter species is deferred 
to a forthcoming article. 


4 CONCLUSIONS 

We have identified the multiple components near [Al VI] as due to 
the HFS splitting of 27 Al. Theory agrees with the observed mag¬ 
netic dipole HFS coupling constants within the uncertainties, giv¬ 
ing support for the use of theoretical constants in the modelling of 
ionic lines profiles. 
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Table 3. HFS fits to the observations 


date, 

noise a 

x 2 /^ 

Ao 

II 

27 Al 

Al,j= 2 

Bl,j= 1 

Bl,j= 2 

fwhmi.AAi 

R 2 , FWHM2 

2003 




fim 

MHz 

MHz 

MHz 

MHz 

kms-\ 10' 

5 fim 

— ,km s -1 

7 May“ 

4.4 

10 -2 

1.02 

3.658715122 

—0 1+ 162 
ui -153 

3489.llg® 

- 0 . 0 !}]} 

6-912 

18-3!;;°,- 

-15.1 

0.56,36.7121 

26 July“ 

4.5 

10 -2 

0.98 

3.6591121 22 

0 o+ 243 

u -231 

3499.ll} 27 

o.ol 2 } 

10.51} 

23.912:2,- 

-5.4 

0.20,65.21“® 

27 July“ 

7.0 

10 -2 

0.87 

3.6593411^ 

417.71 J 2 1 

3332.9l®g 

-54.61“ 

16.912 

99 q+ 0.7 

z,^.o_ 0 7, 

-30.3 

0.26,11.6l2;2 

30 July 6 

6.3 

10 -2 

0.80 

3.6592471}!} 

258.21“° 

3390.2l} 24 

— 19.71}® 

8.71} 

15 9~*~ 0,7 - 
io.y_ 0 8 j 

-18.8 

0.13,8.812;® 

30 July“ 

4.1 

10 -2 

1.08 

3.6591251}} 

276.9lf 0g 

3399. 2+H 

-29.81® 

9.312 

13 1+ 0 - 8 _ 
AO * A —0.8’ 

-9.9 

0.82,31.4l}Jl 

31 July“ 

5.3 

10" 2 

0.68 

3.659365158 

— 117.0!^ 

3575.91^4 

-17.2l}° 

8.812 

17 1+ 0 - 6 _ 
1 • —0.7’ 

-30.2 

0.11,26.7121 



average 


235.5 ±62.8 

3410.1 ±32.4 

-28.0 ±5.3 

10.7 ± 1.0 





“noise used to assess the significance of the fits, normalised to the peak flux density (i.e. S/N = 1/noise). 

6 6-row spectrum centred on the peak of emission. 

“collapsed slit (optimal S/N). 

^Uncertainties on A 0 refer to the last decimal places. A 0 does not match the fine-structure centroid (see text and Ea.171. 
“Uncertainties on all quantities refer to the usual 68.3 % confidence interval (i.e., 1 cr for 1 parameter). 



Figure 4. Bias due to the correlation of the free-parameters for the HFS 
constants Aj = 2 ( 2 ‘ Al) (y-axis) and Aj=i ( 27 A1) (x-axis), from the col¬ 
lapsed spectrum of July 30th. 

We provide the first measurements of electric quadrupole hy- 
perfine coupling constants for any atomic transition in any astro- 
physical object. We discuss the spectroscopic importance of the 
quadrupole terms. The inclusion of the quadrupole terms improves 
the measurement of the magnetic dipole constants, which are oth- 
erwised affected by a statistical bias. 

As an application we have set a 3 a upper limit on the alu¬ 
minum isotopic ratio, 26 Al/ 2 ' Al < 1/33. This is the most stringent 
upper limit on the relative 26 Al abundance in any astrophysical ob¬ 
ject to date. 

However, the accuracy of our measurement is short of quan¬ 
tifying 26 Al production in AGB stars. The expected isotopic ratio 
at the tip of the AGB is at most 1/37, from the ratio of the_^Al 
and 2 'A1 yields in the 6 Mg models of iForestini & Charbonnel! 
C l997i ). The progenitor mass of NGC 6302 is about 5-6 Mg from 
the data summarised in ICasassus et alJ C 200(t) . But the predicted 
26 Al/ 2 ' Al ratio lies below our 3 a upper limit. We can only dis¬ 
card Ri SO = 1 /37 at 2.5 a. Doubling our integration on [Al VI] in 
NGC 6302 would allow a firm test on the theoretical predictions. 


Table 4. The fine structure centroid and limits on the Al isotope ratio. 


date, 2003 

100 x R\ 

100 x R 4 

-''fs 

7 May 3 

1.5±2.5 

4.8± 2.0 

3.659273(35) 

26 July 3 

< 0.67 

-0.4± 2.0 

3.659427(19) 

27 July 3 

1.8±2.6 

-1.5± 2.3 

3.659415(22) 

30 July 2 

1.6±2.5 

2.9± 2.3 

3.659392(22) 

30 July 3 

< 0.88 

-0.9± 1.7 

3.659410(10) 

31 July 3 

< 0.72 

-1.2± 2.2 

3.659427(40) 

average 

- 

0.6± 0.8 

3.659405(7) 


1 formal 1-cr upper limit (see text for accurate upper limits). 

2 , and 3 same as Tabled 

4 best fit isotopic ratio, using fixed HFS constants 
5 rest wavelength in (im and in air. 

To establish useful constraints on the 26 Al production by AGB 
stars we must deepen our observations of NGC6302, and extend the 
analysis to other targets. Only PNe and symbiotic stars have mod¬ 
erate expansion velocities and photoionised coronal-line regions, 
offering narrow emission line profiles in high-excitation species, 
which are otherwise thermally broadened in the Sun. 
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